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ABSTRACT 

We present A^-band (2.12 /im) imaging observations of the SGR 1806-20 field taken 
during its very active phase in mid 2004, which reveal brightening of sources within 
the Chandra X-ray error circle when compared with earlier images obtained in 2002. 

One source brightened by more than a factor of 2, and so we consider this to be the 
probable infrared counterpart for SGR 1806-20. The other two sources are located in 
close proximity to the probable counterpart and show marginal brightening, which may 
suggest that the high-energy photons emitted from the SGR during its active phase 
have induced dust sublimation or brightening of the unresolved background around the 
SGR. 

Subject headings: gamma rays: bursts — stars: neutron — X-rays: stars — pulsars: 
individual (SGR 1806-20) 


1. INTRODUCTION 

The soft gamma-ray repeaters (SGRs) are now commonly believed to be magnetars, isolated 
young neutron stars with ultrahigh magnetic fields of B > 10 14 G (e.g., Duncan & Thompson 
1992). Detections of long period X-ray pulsations and their rapid spin-down rates (P = 5-8 s and 
P ~ 10 -11 ss^ 1 , respectively; e.g., Kouveliotou et al. 1998, 1999; Hurley et al. 1999), cyclotron 
resonance features in the burst X-ray spectrum (Ibrahim et al. 2002), and burst energetics (e.g., 
Paczyriski 1992; Thompson & Duncan 1995) all strongly support the magnetar model. 

Another small class of objects, the anomalous X-ray pulsars (AXPs: Mereghetti & Stella 1995), 
are also presumed to be magnetars, owing to their similarities with SGRs both in spin period and 
deceleration rate. In addition, the discovery of several SGR-like X-ray bursts from AXPs (Gavriil 


1 Based on data collected at the Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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et al. 2002; Kaspi et al. 2003) has increased the evidence for a connection between SGRs and 
AXPs. To date, possible optical-infrared (optical-IR) counterparts with unusual colors have been 
found for five out of six AXPs. Four show variability in some way correlated with changes in 
X-ray flux (4U 0142+61: Kern & Martin 2002; Hulleman et al. 2004; IE 2259+586: Kaspi et al. 
2003; IE 1048.1-5937: Israel et al. 2002; and XTE J1810-19704: Rea et al. 2004). In the case of 
AXP IE 2259+586, the iv s -band flux increased by a factor of ~ 3 by the third day after an X-ray 
burst. Thus, variability is thought to be a common characteristic of the optical-IR counterparts for 
AXPs. In contrast, none of the five SGRs have convincing counterparts detected with optical-IR 
variability, despite considerable observational effort. The detection of optical-IR counterparts for 
SGRs is essential in order to investigate the connection between SGRs and AXPs. If optical-IR 
counterparts for SGRs show a relation between variability and burst activity similar to the AXPs, 
it will enable much easier detection of SGR counterparts after outbursts. 

SGR 1806-20 gradually entered an active phase in late 2003 (e.g., Hurley et al. 2003). It showed 
recurrent intense bursts in mid 2004, and a series of many bursts in late 2004. Finally, a giant flare, 
the first observed from SGR 1806-20, occurred on 2004 December 27 (e.g., Borkowski et al. 2004; 
Hurley et al. 2005), and left a fading radio afterglow with a precise localization (e.g., Cameron et 
al. 2005; Gaensler et al. 2005). A series of infrared observations were performed with the Subaru 
8.2-rn telescope during the active phase of SGR 1806-20 in mid - late 2004, several months prior to 
the giant flare. Here we report the discovery of the probable infrared counterpart for SGR 1806-20 
within the Chandra X-ray error circle. We also discuss the influence of the high-energy photons 
emitted by the SGR burst upon the surroundings. 


2. OBSERVATIONS, DATA ANALYSIS, AND RESULTS 
2.1. Infrared Data and Photometry 

Our first epoch observations were obtained on 2002 May 25 in J (1.25 /im), H (1.63 /irn), 
and K' (2.12 /im) bands using the infrared camera and spectrograph (IRCS: Kobayashi et al. 
2000) on the Subaru 8.2-rn telescope (lye et al. 2004) under photometric sky conditions. A plate 
scale of 07058 pixel -1 was chosen so as to match the average seeing conditions (FWHM = 0730 
in iC-band) on the night of the observation. The field of view (FOV) was approximately 1' x 1'. 
Images were taken by using nine-point dithering of three coadds with 30-s exposures each, and the 
resulting total integration time was 810 s for all three bands. A UKIRT faint standard star FS 
148 (Casali & Hawarden 1992) was also observed using the same configuration. We carried out 
a second epoch of observations in the K' band on 2004 August 8, during the active phase, with 
the same configurations as for the first epoch, to allow a precise comparison of the images from 
the two epochs. The sky conditions were again excellent, and the seeing size was 0732. A total 
integration time of 1260 s was achieved with a nine- and a five-point dithering of 90-s exposures at 
each position. We also performed another observation in the K' band with adaptive optics (AO: 
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Takami et al. 2004) on the same night using a 13.7-th magnitude star located 30” north of the 
SGR as the natural reference star. The plate scale was set to 0 / . , 023 pixel -1 (FOV was 23" square), 
and the stellar image size improved to 0720. Two sets of nine-point dithering observations of 90-s 
exposures at each position were obtained. Fourteen frames with good seeing were selected, and the 
resulting total integration time was 1260 s. The 3rd and the 4th epoch observations were carried 
out using AO on 2004 September 3 and 23, respectively, also during the active phase. The best 
quality frames with three coadds of 30-s exposures were selected to create final images for those 
epochs. The resulting total integration time and the stellar image size were 1710 s and 0722 for the 
3rd epoch, and 1980 s and 0738 for the 4th epoch. 

The standard imaging calibration procedures were applied for each epoch data using the IRAF 
data analysis software. 2 There is some latency in the ALADDIN array when bright objects are 
focused on it; therefore, the locations of bright objects in the previous frame were masked before 
stacking selected images. Then, a photometry of the field stars in the vicinity of the SGR X-ray 
position was carried out with the DAOPHOT/ALLSTAR package (Stetson 1987) in IRAF. The 
radius of the point spread function (PSF) was set to twice the stellar FWHM at each epoch, and 
the aperture radius for the photometry was set equal to the FWHM. We applied the photometric 
standard star FS 148 to our first epoch data and measured a field star C (see Figure 3 in Eikenberry 
et al. 2001, and Figure 2D in this paper) as a secondary standard. The Id-band magnitude of star 
C was I\' = 16.48 ± 0.05, which is consistent with previous work by Wachter et al. (2004) but 
slightly different to that of Eikenberry et al. (2001). For the other epoch data, photometry was 
conducted relative to star C, based on the assumption that the star is not variable. The 3 -a limiting 
magnitudes were 21.5, 21.4, 21.9, 22.0, and 21.0 for the 1st epoch, 2nd epoch without AO, 2nd 
epoch with AO, 3rd epoch, and 4th epoch data, respectively. 


2.2. Optical Data and Astrometry 

An unfiltered optical image of a 10-s exposure, covering a FOV of 6 7 diameter, was taken with 
the faint optical camera and spectrograph (FOCAS: Kashikawa et al. 2002) on the Subaru telescope 
on 2002 June 18. As the seeing was 079 and the limiting magnitude was sufficiently deep, we used 
this image as a reference to obtain precise astrometry for our IR data. First, the positions of 22 
unsaturated, stellar-like sources in the optical image were measured and compared with the USNO- 
B1.0 catalog (Monet et al. 2003). The optical image was then compared with the J-band image 
using 20 stars common to both images. Finally, the J-band image was compared with the JV-band 
image using another 16 stars. Taking the positional errors of the individual stars from the USNO- 
B1.0 catalog into account, the astrometric errors incurred in these procedures were 0706, 0701, and 
less than 0701, respectively, in each coordinate. An empirical value for the size of local offsets in the 


2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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southern part of the USNO-Bl.O catalog compared to the International Celestial Reference Frame 
(ICRF) is O'/136 (de Silva Neto et al. 2005). The resulting final astrometric uncertainties in our 
R'-band images are tf'15. Thus, the 99% uncertainties for the Chandra X-ray error radius (O'/S; 
Kaplan et al. 2002) and the VLA radio afterglow error radius (0 / . / 1; Cameron et al. 2005) on our 
infrared images are 0 / /TT and 0 / /41, respectively. The positions of the field stars, along with the 2nd 
epoch photometry, are summarized in Table 1. 


2.3. Field Around SGR 1806—20 

The SGR 1806-20 field along with its putative associated massive cluster is shown in Figure 
1A. A luminous blue variable star LBV 1806-20 is seen 12" east of the center. Bright cluster stars 
are located mainly to the north to north-west of the SGR position. As mentioned by Corbel & 
Eikenberry (2004), the color of these stars indicates that they suffer large extinction {Ay ~ 30). 
A region of even higher reddening can be observed about 25" to the north-west of the center. 
The close-up (Figure IB) shows a lot of faint and redder, consequently probable cluster member, 
stars. There are two remarkable moniliform stellar arcs 2'.'l north and 3'.'6 south-west of the SGR. 
Unfortunately, ghost images of the bright star, produced by the compensator inside the IRCS optics, 
overlap with part of the stellar arc to the south-west. It is intriguing that the center of curvature 
of both arcs is located close to the SGR position. 

Multi-epoch /i'-band images are presented in Figure 2. Figure 2D shows the Chandra and 
VLA positions with 99% confidence error circles (see 2.2) along with the labels of individual field 
stars. Stars A, B, and C in Eikenberry et al. (2001) correspond to our stars Al to A4, B1 to B3, 
and C, respectively. Better spatial resolution and a deeper detection limit enabled us to spot the 
multiplicity of the stars A and B. A remarkable difference can be seen between Figure 2A and 2B: 
star B3 can hardly be seen in the 1st epoch but is bright in the 2nd epoch. There is no other major 
difference between the two images. All good quality images taken with AO are accumulated to 
create a high-resolution image (Figure 2C; FWHM = O'.'18). 

The light curve of each field star is presented in Figure 3. Most of the stars are constant within 
the error bars, which implies that star C, used in calibration, was constant during the period covered 
by our observations. B3 shows brightening by more than a factor of 2, whereas A3 shows a small 
amount of fading. B1 and B2 show slight brightening that just exceeds the error bars. 


3. DISCUSSIONS AND CONCLUSIONS 
3.1. Detection of a Variable IR Counterpart 


The position of SGR 1806-20 is well localized in the X-ray; however, the held of SGR 1806- 
20 is very crowded at IR wavelengths and the probability of a chance coincidence with unrelated 
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IR sources is relatively high. Therefore, variability is a key means of distinguishing the real IR 
counterpart. The IR source B3 is the only source to show brightening in the K' band clearly during 
the active phase. The chance probability of finding a variable star is generally low, i.e., 4.1 x ICR 3 
is derived from the Optical Gravitational Lensing Experiment (OGLE-II) dataset (Zebruri et al. 
2001), and is much lower if we select highly variable stars like B3. 

Recently, the persistent hard X-ray emission (<200 keV) from SGR 1806-20 was presented, 
and its spectral shape was different from that of the bright bursts (Molkov et al. 2004). It is also to 
be noted that the persistent hard X-ray emission correlates in intensity and spectral hardness with 
the level of bursting activity (Mereghetti et al. 2004). The flux in the 20-100 keV energy band was 
almost doubled during the active phase in 2004 September to October compared with the data in 
2003. On the other hands, our first epoch data was taken in 2002, and the IR flux for B3 increased 
by more than a factor of 2 in our third epoch in 2004. It is likely that the IR flux would increase 
with the persistent X-ray flux. 

We should also note here that two short bursts were detected on 2002 August 25 by Konus- 
Wind and Ulysses (Hurley et al. 2002), just 1.5 and 9 hours prior to our iG-band observations. Our 
second epoch observations were made 11.4 and 15.9 days after the short bursts on July 28 and July 
23 with fluences of 4.4 x 10 -6 erg cm -2 and 7.5 x 1CU 6 erg cnr -2 , respectively (Golenetskii et al. 
2004a). Our third epoch data were collected 5.8 and 8.3 days after two intense intermediate bursts 
on August 28 (4.0 x 10 -5 erg cm -2 ; Golenetskii et al. 2004c) and August 25 (2.1 x 10 -5 erg cm -2 ; 
Golenetskii et al. 2004b), respectively. Some intense SGR bursts are followed by an X-ray tail or 
afterglow that decays with time (e.g., Ibrahim et al. 2001; Lenters et al. 2003). Assuming that 
the SGR has an IR afterglow associated with X-ray burst activity, by analogy with AXPs (e.g., 
AXP IE 2259+586: Kaspi et al. 2003; Tam et al. 2004; and XTE J1810-197: Rea et al. 2004), SGR 
1806-20’s IR counterpart would be expected to be brighter at the third epoch than at the second, 
because the elapsed time since the last burst was shorter and the fluence of the burst was larger. 
The light curve of B3 between the 2nd and 3rd epochs is also consistent with this model. On the 
other hand, the elapsed time period since the last burst for the 1st epoch data is much shorter than 
for the 2nd epoch data, but B3 is only marginally detected in the 1st epoch. This implies either 
that 1) time lags exist between X-ray and IR flux enhancements, or 2) multiple burst activities are 
in some way accumulated in the IR. 

An alternative interpretation of the brightening behavior of B3 is that B3 could have a large 
proper motion (0 / /l yr -1 ) and so was very close to B1 in the 1st epoch. In this scenario the 
difficulties in accurately separating B1 and B3 may have caused the apparent brightening of B1 
and fading of B3 at the 1st epoch. Future high-resolution imaging observations are necessary in 
order to measure precisely the proper motion of these sources. 
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3.2. Brightening of Nearby Stars 

The active phase of SGR 1806-20 has been ongoing for several months, and so the inter-stellar 
medium (ISM) surrounding the SGR may be affected by the presence of high energy photons. In 
the case of gamma-ray bursts, for example, a large amount of dust sublimation is expected by 
the optical-UV prompt emission (Waxrnan & Draine 2000) and X-ray radiation (Fruchter, Krolik, 
& Rhoads 2001). It is suspected that SGR 1806-20 is associated with a massive molecular cloud 
(Corbel & Eikenberry 2004), and suffers a large extinction (Ay > 30) inferred from its X-ray column 
density of Njj ~ 6 x 10 22 cm -2 (Sonobe et al. 1994; Mereghetti et al. 2000) and the relation of 
Ay I Nh = 5.6 x 10 -22 (Predehl & Schmitt 1995). As the burst activity has been ongoing for a 
considerable period, the dust in a spherical region centered on the SGR is expected to have been 
sublimated, and hence the brightness of the background stars close to the line of sight of SGR 1806- 
20 might be expected to brighten. In fact, we marginally detected subtle brightening of the nearby 
stars B1 and B2. However, the brightening of the nearby stars could also be interpreted as the 
brightening of unresolved background (nebulosity) owing to the heating of the dust surrounding the 
SGR. We can discriminate between these two mechanisms by measuring the variability in different 
wavelengths. In the first scenario, the color of the stars may change, i.e. they will be brighter 
at shorter wavelengths, and the scale of the variability follows the extinction law. In the latter 
case, the variability may be larger at longer wavelengths. The distances of the stars B1 and B2 
from the line of sight toward B3 are fy/23 and (//30, which correspond to 20.3 and 26.0 light-days, 
respectively, at a distance of 15.1 kpc (Corbel & Eikenberry 2004). Therefore, the burst photons 
could clearly have reached and affected the ISM along the line of sight towards B1 and B2 when the 
observations were obtained. Further multicolor monitoring is required to confirm these conclusions. 


We would like to thank all staff members of the Subaru Telescope for their observational 
and data analysis support, and the FOCAS team members for the use of the data taken in their 
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Fig. 1.— Wide-field (55 w x 55 w ) J, H, and /i'-band three-color composite image of the SGR 1806- 
20 field taken on 2002 May 25 (left). Blue, green, and red correspond to the J, H, and K' bands, 
respectively. North is Tip and East is left. The elongated sources at north and west from bright 
stars are ghost images caused by the beamsplitter and compensator, respectively. A close-up of the 
central region (gray square with 10 // x 10 // ) is shown in the right panel. A gray circle represents 
the Chandra X-ray error circle (Kaplan et al. 2002) with a radius of 0 / /77 (see 2.2). A compensator 
ghost is seen to the south-west of center. 
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Fig. 2.— (A) Ji'-band image taken on 2002 May 25 (1st epoch) without using AO. (B) iF'-band 
image taken on 2004 August 8 (2nd epoch) without AO. (C) A''-band image taken with AO. Only 
images with smaller FWHM in the 2nd and the 3rd epoch are combined. The final stellar size is 
O'.'18 in FWHM. (D) Same image as (C) is shown with a different intensity scale, showing the 0'.'77 
radius error circle around the Chandra X-ray position of a = 18 /l 08 m 39 s .32 and 6 = — 20°24'39'.'5 
(J2000.0), the O'.'41 radius error circle around the VLA radio afterglow position of a = 18 ft 08 m 39 s .34 
and 6 = —20°24'39'.'7 (J2000.0), and the star labels in Table 1. 
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Fig. 3.— K '-band light curves of the SGR 1806-20 held stars are presented. Solid lines are used 
for the stars A1 - A4, B1 - B3, and C, and dashed lines for the stars FI - F7. Owing to the larger 
plate scale, we were unable to separate star A4 from A1 in the O'/OSS pixel-scale images (without 
AO). Therefore, the flux from A4 contaminates the flux of A1 in the 1st epoch, and hence the 
magnitude of A1 may be overestimated by about 0.1 mag. B3, F7 and A4 were not measured in 
the 4th epoch data owing to the relatively poor seeing conditions. The 2nd epoch plots are based 
on the observation with AO. 
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Table 1. Astrometry and Photometry of the Point Sources Near SGR 1806-20 


Label RA(J2000) DEC(J2000) K' magnitude 


Al 

18 

08 

39.386 

-20 

24 

39.02 

17.63 

± 

0.05 

A2 

18 

08 

39.369 

-20 

24 

39.43 

17.91 

± 

0.05 

A3 

18 

08 

39.358 

-20 

24 

39.01 

19.27 

± 

0.07 

A4 

18 

08 

39.383 

-20 

24 

39.18 

19.88 

± 

0.08 

B1 

18 

08 

39.340 

-20 

24 

39.77 

19.39 

± 

0.07 

B2 

18 

08 

39.343 

-20 

24 

40.14 

19.36 

± 

0.07 

B3 

18 

08 

39.329 

-20 

24 

39.94 

21.00 

± 

0.36 

C 

18 

08 

39.319 

-20 

24 

40.81 

16.48 

± 

0.05 

FI 

18 

08 

39.382 

-20 

24 

40.34 

19.72 

± 

0.08 

F2 

18 

08 

39.290 

-20 

24 

40.31 

19.32 

± 

0.07 

F3 

18 

08 

39.277 

-20 

24 

39.26 

20.00 

± 

0.12 

F4 

18 

08 

39.257 

-20 

24 

39.24 

19.52 

± 

0.07 

F5 

18 

08 

39.214 

-20 

24 

40.27 

19.19 

± 

0.06 

F6 

18 

08 

39.229 

-20 

24 

40.85 

19.45 

± 

0.06 

F7 

18 

08 

39.223 

-20 

24 

38.63 

20.58 

± 

0.13 


Note. — Measurements are based on the dataset taken with AO on 2004 August 8. Units of 
right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, 
and arcseconds. Astrometric uncertainties are 0 , . / 15 in each coordinate. 
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